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ON THE METHOD OF CONTINUED IDENTITY.* 

By Mb. Chas. H. Kummell, Washington, D. C. 

This apparently new method of expressing a quantity is based on the fol- 
lowing considerations : — 

Suppose a quantity x is approximately known by its functional properties, 
so that we can place 

where l/g is the approximate value of x and x^/q its correction. Again, 
assume 
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Suppose any one x, for instance x^, =0, or sufficiently nearly so; then 
we have 



X = 



1 + «1 



1 + « 



V (1 + <^2) 



= « , a perfect identity ; 

and it is important to notice that this identity holds, whatever the preceding 
values of a^g and x^ may have been, if only ajg = or is sufficiently small to be 
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neglected. We may thus continue the identity to any number of steps before 
perfect or sufficiently perfect identity is reached. Since 

1 -h «i 1 • i. 1 
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the condition for perfect identity of (2) is that the //s must have co for their 
limit. Take, for instance, the binomial series which may be written 

1 + . ■ . ; (4) 



(5) 
(5x) 

(5,) 

The condition for perfect identity is then 

g'^ = 00 ; .-. rn — n -{- 1 =^ 0. 

It follows, then, that the binomial series can give perfect identity only if m is 
an integer, and that this is attained by m + 1 terms. If m is not an integer, 
perfect identity is impossible, since ^„ can never be infinite ; and this is also 
the case if the condition of convergence a; < 1 is fulfilled. Since, then, 



{i + xr=^i + ''^x 
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comparing which with (2), we have 
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the terms at infinity form a geometric series with the exponent — x, which 
therefore has a sum. It follows, then, that the form of continued identity is 
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still convergent if the denominators q tend to a limit > 1. This condition 
holding, we can give an approximate value of the remainder if we stop at any 
term ; for if y„ is sufficiently large, the next term, though small, will exceed the 
sum of all following terms ; and since </„^.j must be at least = <?„, we have 

(?) + l)th term <— . — . — . . . — = approximate remainder. (6) 

<I 'll 'Ji In <ln 

Transforming any equation for one unknown quantity x by relations (1) 
into an equation for x^ , x^ , . . ., we have a method of approximation which 
is convenient and, in general, rapidly convergent. Suppose we have to solve 

=/(»•). (7) 

Let 

X = a;,, + a?i , (8) 

where x^ is the integral and x^ the fractional part of the root. By Taylor's 
theorem we have 

=/(^o) +/>o) ^t 4- ^^ «•! ' + ^^ x,' + . . . . (9) 

Let 



(1,) 

fly, \ 

where q is the integral part of — ^)^[ (preferably exceeding it). To transform 

(9) to an equation in a- 2 we have to decide first to which order of precision we 
wish to satisfy the equation. To fix ideas, let us suppose the solution to be 
exact to the 3d order ; we have then to solve a cubic equation, which we shall 
write thus :■ 

■=: u^'" + a/'ajj + aj'j?!^ + a-^x^^. (9i) 

Transforming this to one having its roots q^ times as great, we have 

then transforming this to one having its roots diminished by 1, we have 
= «i +"iVi +«i"2'i^ +«i"'!?* 
= a^"' + a^"x^ +a^x^^ +a^x^^. (Bj) 

Let ^, = ^^±^, (1.) 
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where 172 = integral part of ^y,-; (preferably exceeding it). Transforming to 

12^2 — 1 *s before, we have 

= «3"' + «3"«3 +«3'«3^ +«3«3*- 

This process is continued as far as we please, in case the cubic (9 J is 
proposed for solution. If, however, (9^) coincides with the complete equa- 
tion (9) to the 3d order only, it will be useless to make more than three trans- 
formations. A convenient form for making these transformations I shall 
exhibit in the solution of the cubic 

= 3 — 5a; 4 x\ 

This equation has three real roots, one between and 1, another between 
1 and 2, and the third between — ■ 2 and — 3. Here f is an approximate value 
for the smallest root ; hence /^ = 2 > f, and a more exact value of the root is 
found by the following algorithm : 



0=r:+ 3 — 5x + Ox^ +x^ a 

+ 6—10 hi 

+ 12—20 +1 +1 

+ 24—19 +2 +1 

0=+ 5 — ITa-'i + 3.z'i2 +«,» x^^ — 



. 1 +a?i 



+ 20—68 +12 +1 

+ 80—272 +13 +1 

+ 320 — 259 +14 +1 

0=+ 61 — 245a;3 + ISa!^^ + x^^ a 

+ 244 — 980 +60 +1 

+ 976 — 3920 +61 +1 

+ 3904 — 3859 +62 +1 

0=+ 45 — 3797a;3 +63a;32 + x^^ X3=^'^'"* 



4 " 



84 

The root so far determined is correct within less than 7, — r — a — rm! = stvpjttt- 

2 X 4 X 4 X 84^ 86016 

We have, then, neglecting x^, 

x^ = 0.011905, 
x^ = 0.252976, 
«! =^ 0.313244, 
X =0.656622. 
The correct value is x = 0.656620. 
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To determine another root, for instance that lying between — 2 and — 3, 
we may first increase the roots by 2 by assuming x = y — 2 ; then y will lie 
between and — 1, and the more exact value is found thus : 

0=+ 5 4- ly - 6y^ +y^ y =^^ 

— 10 — 14 +12 +1 
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= _ 824 +1647282/3 -822^^2 -^y^^ y^^lAilA 
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This determination must be correct within less than 



2 X 55 X 2002 4400000' 



and we have 



j/j = + 0.005, 
j/i = — 0.0182727, 
y =— 0.4908636, 
X = — 2.4908636. 

The correct value is a? = — 2.490863S. 

We notice the greater convergence in the computation of this root. This 
is caused by its being more different from the other roots. 

Should two or more roots be very close together, the process becomes 
very slowly convergent. This is, indeed, a difficulty common to all numerical 
methods of solution of higher equations. It may be in most cases sufficiently 
overcome by placing the nearly equal roots between and 1, and then deter- 
mine their reciprocals. For a test of this suggestion, I solve the equation, 

= 7 — 7a; + a!^ 
which has two roots between + 1 and -|- 2. Here the direct process gives 



1.3568948. 



If, , ir, , If, 1 



l + 4[l + 3tl + 3ll-? 



1 + B 
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If we assume x ^y -\- 1, then 

= 1 - % + 32/2 + 2/^ 
and placing y ^ s-^, we have 

= 1 +dz — ig^ + 3^. 

This equation has one positive root between + 1 and + 2, and one be- 
tween + 2 and + 3. To determine the larger, which corresponds to the root 
already given, we assume z = v -\- Z, and have 

= + 1 + Qv + 5w« + v^ 

_ (5 _ 36 — 30 +1 

+ 36 + 216 — 29 +1 

— 216 +187 — 28 +1 
= — 29 + 159«;i — 27^1 2 +v-,^ 

— 145 + 795 — 135 + 1 

— 725 +3975 —134 +1 

— 3625 + 3841 — 133 + 1 
0=+ 216 +37081)2 — 132«22 _|_^^a 
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This should give the root within less than ^ = =-=^ = -.._„ , and we have 

X o X 1' oo70 

v^ = — 0.05882, 
wj = + 0.18824, 
V ^— 0.19804, 
3 = + 2.80196, 
y = + 0.35689, 
« = + 1.35689. 

The exact value is « = + 1.356896. 

If we apply the algorithm to the equation 

= 8 — 7a! + «% 

which differs from the above equation by 1 in the numerical term, the two 
roots corresponding to those before nearly equal will be imaginary, and we 
have 
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We have, therefore, the periodic form 
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which obviously has no value. 

The application of the method to the general quadratic 







is of great interest. Assuming x - 

tion in it'i 

= 1 + (2 



-\- a'x + x 
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L , we have the transformed equa- 
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Assuming 



we have the equation in a;^ 

= 1 + [2 — {a'^a"-' — 2)2] x^ + a;^^ 
We have, therefore, 
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or if we place 
then will 

We have 
therefore, 



a 'a 



1 + 



q = — aa ~, 

q^ = — a'^a"~' — 2, 

?2=?i' — 2, 
?3=?8' — 2, 



oe = 






\T " ~2 « 
Placing 

and dividing by ^ , 



1+-7, 



1 + 



?.'-2 



1 + 



1+. 



]a+3=--l + 



2[^ 4s-' + 2l 



4a" 



1 + 



16s- 



163- 



1+. 



(10) 



(11) 



(12) 
(13) 

(14) 



Now this form, which is essentiallj' Kobert's rule for extracting a square root,* 
has the remarkable property of being convergient for any value of z, while its 
value by the binomial series 
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* Ooiiunuuieated to the Mathematical Section of the Philosophical Society of "Washington by 
Dr. A. Martin, and published in the Mathematical Magazine, Vol. II, p. 13, as problem 132. 
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ceases to be convergent if s > 1 ; for the denominators q, q^, q^, . . . can 
never be < 2, whatever the value of s. If, however, s < 1, these denomi- 
nators tend to infinity, which will bring about perfect identity, while the form 
is only convergent if s > 1. 
If we assume 



f{x)=. 



1 + 



^i*" 



■?i L 



2«~^+ 2'2'«"'+2'2 



?3«~' + ?3'«"' + ^s"*"' + ?3 



1 + . . . , (15) 



we have a form which becomes perfectly identical if « < 1, and is convergent 
only if a; > 1, provided the quantities q, q^', q^", q^'", . . . tend to values > 1. 
We may assume identically 
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and by McLaurin's theorem we have 

fM =/i(0) +/i'(0) . X +'^^x' + 



hence, if/, (0)=0, 
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We have further, if /^(O) =/2'(0) == 0, 
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and so on. 

If the proposed function f{x) vanishes for a; ^ by virtue of a factor 
x" (where p may be either whole or fractional), we apply the above process 
to the function 

x~^f(x), denoting it f(x), 



or to 



x-^^'f(x). 



„-p + 2 






If the function is par, or /(«) =/(—«), its development contains only 
even powers of x, and placing a;^ = z we have the required form. If the 
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function is impar, or f{x) = — /(—«), dividing by x and placing x^ = s -we 
have again the required form. 

We may also arbitrarily modify the complete form (15), if, for instance, 
we do not determine all the ^''s at the nth term. This, however, changes noth- 
ing essential in the rule in equations (18) and (19). At any step there will be 
as many equations (18) as we choose ; then formulae (19) follow where (18) have 
been interrupted. To explain this modification, let us suppose we only want 
two ^'s in the fourth term, then we have 

0=/4 (0) = -l+^3-^-^ •••!73=j4|oy. (I83) 

o y'/n^ yV^'W , '/3"'(0) > 3! /3'^(0) ,.0 « 

=/4 (0) = ?3 -^'-jf^ + ?3 "^'gH ' ■'■ ?3 = "~77X0) ^8 4!~ ' (^^3 ) 
/4"(0)_^ /3^(0) , ,/3'^(0) 

-21 ^8 5! ^ ^3 41 ' 

A"'{0) _ ^ /3'"(0) I ^ ^/3^(0) ng) 

~3! ^' 6! ^ ^' "TT ' ^ ' 



It happens in some cases, much to the advantage of the form (15), that one or 
more of equations (19) vanish. In that case a corresponding number of terms 
in (16) also disappear. 

For the first application take the function 
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We have then 

o = !Zi/i'(0)-i; .■.q, = i, 
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Since here ■ ^n) - vanishes, the third term then disappears, and we have 
o ! 
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^ = 700^"-^' •••?4--3-, 
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etc., etc. 

"I ft 

Since here ^1' = ^, q^" =^, ?/ —2673' 
a value > 1, we cannot assert that the form 



1 7700 J^ 1400 3^ 2450 
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1 1400 , 1 2450 
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do not appear to tend to 



1 + 1- 
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1 + 
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12 + 12* + -"«= 
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700 ,1400 ,7700 3 , 1400 3 , 2450 
"3~"^~3""'^'2l^'' "•"'27"" "'"2673' 



(20) 



is convergent for every vahie of x. It may, however, be used for values of 
a > 1, to a considerable extent, as long as the denominators q, ^i«"^ + q^, 
q2^~^-\- qi^~^ + ^s") • • • form an ascending series. If a; < 1, the form is 
perfectly identical, and the first two terms will be correct to at least six deci- 
mal places. 

Placing X = , we have 



l{l-\.x)=l 



1 + 



= -^1 



98 



KUMMELL. ON THE METHOD OF CONTINUED IDENTITY. 



and by (20), 

1(1 -{-«)= — Z (1—2) 
+ 



1-2^ 



1 + 



12 — 123 + 1^2 



Since z 



700_1400 ,7700 ^ 1400 3 2450 4 

-3- -3- 3 -t- 27 ^ 27 ^ "*" 2673 ^ 

- < 1, this form is always perfectly identical. 



(21) 



This example is a fair exhibit of the advantages and defects of the 
method. Although the required coefficients g are determined by a regular 
algorithmic process, the resulting form is by no means of such a form that 
it might be continued by induction, except in case of Robert's rule (14) and 
the formula for the root of a quadratic (10). 



